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Abstract. For the Toeplitz quantization of complex- valued functions on a 2n-dimensional 
torus we prove that the expected number of eigenvalues of small random perturbations of 
a quantized observable satisfies a natural Weyl law (1.3 1. In numerical experiments the 
same Weyl law also holds for "false" eigenvalues created by pseudospectral effects. 



1. Introduction and statement of the result 



In a series of recent papers Hager-Sjostrand [15] . Sjostrand [T7], and Bordeaux Montrieux- 
Sjostrand [3J established almost sure Weyl asymptotics for small random perturbations of 
non-self- adjoint elliptic operators in semiclassical and high energy regimes. The purpose of 
this article is to present a related simpler result in a simpler setting of Toeplitz quantization. 
Our approach is also different: we estimate the counting function of eigenvalues using traces 
rather than by studying zeros of determinants. As in [T3] the singular value decomposition 
and some slightly exotic symbol classes play a crucial role. 

Thus we consider a quantization C°°(T 2n ) 3 f i — > /jv G Mjyn(C), where T 2n is a 2n- 
dimensional torus, R 2n /Z 2n , and M N n(C) are N n x N n complex matrices. The general 
procedure will be described in ^2] but if n — 1 and T = § x x§^, then 

/ = f{x) .— f N = diag (f(£/N)) , € = 0,.-. ,JV-1, 

g = g(0^9N = J r * N ^g(g(k/N))^ N , k = 0, ■ ■ ■ , N - 1 , 

where JF/v = (exp(2irikl /N) / \ZN) <k,e<o,N-i, is the discrete Fourier transform. 

Let uj i — > Qn(u) be the gaussian ensemble of complex random N n x iV n matrices - see 
£j3j With this notation in place we can state our result: 

Theorem. Suppose that f G C°°(T 2n ) ; and that ft is a simply connected open set with a 
smooth boundary, dQ, such that for all z a neighbourhood of dVt, 

(1.2) vo\ T 2n({w : \f(w)-z\ <t}) = 0(t K ), 0<t<l, 
with 1 < k < 2. Then for any p>pQ>n + l/2 

(1.3) E w (| Spect/W + iV p g7v(^)) n ft|) = N n vo\ T 2 n (r\Q)) + G(iV ra -' 3 ) , 
for any (3 < (k — + 1). 

l 
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Figure 1. On the left we reproduce [7J Fig. 1.2] with two types of regions 
used for counting added. It represents Spec(/iv+£') where N = 100, f(x, £) = 
cos(27nr) +icos(27r£) (called the "the Scottish flag operator" in [7]), for a 
hundred complex random matrices, E, of norm 10~ 4 . On the right we show 
the counting functions for the two regions, and the corresponding Weyl laws, 
as functions of r. The breakdown of the Weyl law approximation occurs when 
the norm of the resolvent (/# — -z) -1 , \z\ = r, or | Rez| = r, is smaller than 
WE]]- 1 = 10 4 . For tt = {\z\ < r}, r < 1, k = 2 and for Q = {\Rez\ < r}, 
k = 3/2 at four points of dfl (intersection with the boundary of /(T)). For 
r = 1, the corners satisfy (1.2) with k — 1. 



Remark. The theorem applies to more general operators of the form A(N) = + g^/N, 
where g may depend on iV but all its derivatives are bounded as iV — >• oo. 



The main point of the probabilistic Weyl law (1.3) is that for most complex- valued 
functions / the spectrum of fa will not satisfy the Weyl law. Yet, after adding a tiny 
random perturbation, the spectrum will satisfy it in a probabilistic sense. As illustrated 
in Fig|2] a tiny perturbation can change the spectrum dramatically, with the density of the 
resulting eigenvalues asymptotically determined by the original function /. 



Condition (1.2) with < k < 2 appears in the work of Hager-Sjostrand |13j . Its main 



role here is to control the number of small eigenvalues of (/jv — z)*(fa — z), see Proposition 



2.9 and that forces us to restrict to the case k > 1. It is a form of a Lojasiewicz inequality 
and for real analytic / it always holds for some k > 0, as can be deduced from a local 
resolution of singularities - see [U Sect. 4]. Similarly, for / real analytic and such that 
f(T 2n ) C C has a non-empty interior, 



df\ f -i( z)J L ==> ([L2J) holds with k > 1 
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Figure 2. The MATLAB computed spectra of f N for f(x,g) = cos(27ra) + 
icos(27r£). For N = 100 the computations return the correct eigenvalues 
following the Scottish flag pattern. For N = 1000 the actual spectrum of 
still follows the same pattern but the computations return "false eigenvalues" 
which satisfy the same Weyl law as random perturbations. The plots of 
the counting function for a single random matrix are very close the Weyl 
asymptotic even in the case of N = 100 providing support for the conjecture 



For / G C°°(T 2n ) we have 



(1.2) holds with k = 2, 
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and by the Morse-Sard theorem the condition on the left is valid on a complement of a set 
of measure in C. Also, 

\/w e f-\z) {/, f}(w) ± or {/, {/, f}}(w) ± (Q holds with « = |, 



see [T3l Example 12.1]. Here {•, •} is the Poisson bracket on T 2n (see (2.22) below). 

The significance of the Poisson bracket in this context comes from the following fact: 

(1.4) {Re/,Im/}H<0, z = f(w) \\(f N - z)'^ > N*/C p , \f V > , 

and moreover an approximate eigenvector, w^v, causing the growth of the resolvent can 
be microlocalized at w (meaning that for any g vanishing near w, \\gNUN\\p = O^N^ 00 ), 
1 1 un 1 1 e 2 — 1, see ^2J). This is a reinterpretation of a now classical result of Hormander proved 
in the context of solvability of partial differential equations - see [8], [20], and references 
given there. For quantization of T (1.4) was proved in [7J, and for general Berezin-Toeplitz 
quantization of compact symplectic Kahler manifolds, in [5]. 



The relation (1.4) shows that {/, /} 7^ implies the instability of the spectrum under 
small perturbations. In that case the theorem above is most interesting, as shown in 
Figures [l] and |2} However, as stressed in [3], [13], and [TT] , the results on Weyl laws for 
small random perturbations have in themselves nothing to do with spectral instability. For 
normal operators they do not produce new results compared to the standard semiclassical 
Weyl laws: the distribution of eigenvalues is not affected by small perturbations and satisfies 
a Weyl law to start with. 

The numerical experiments suggest that much stronger results then our theorem are true. 
In particular we can formulate the following 



Conjecture. Suppose that (1.2) holds for all z e C with a fixed < k < 2. Define random 
probability measures: 



N r> 

zeSpec (f N +N-PR N (uj)) 



Then, almost surely in u, 



H N {u) — >f*(a n /n\), N — > 00 , 
where o = J2k=i d£k A dx^, € T 2n , is the symplectic form in T 2n . 



The result should also hold for more general ensembles than complex Gaussian random 
matrices. Sjostrand's recent paper [17] suggests that random diagonal matrices would be 
enough to produce the Weyl law-creating perturbations. 

Bordeaux Montrieux [2] pointed out to us that by taking singular /'s, or /'s for which 



derivatives grow fast in iV (corresponding to p = 1 in the S p classes described in £2.1), usual 
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Toeplitz matrices fit in this scheme and that numerical experiments indicate the validity 
of Weyl laws in this case. 

Hager [T2] indicated how the methods of [32] should apply to the case of Berezin- Toeplitz 
quantization but that approach did not suggest any simplifications in the method. In this 
paper we follow the most naive approach which starts with the following false proof of the 
theorem: 



Spec (f N ) n fi| = / tr(f N - zY l dz 

"=" Nn w~ I I (/H - z)- 1 dC{w)dz + o{N n ) 

^ m JdQ Jj2n 

= N n [ (— [ (f(w) - z)- l dz) dC(w)dz + o(N n ) 
J T 2n \2-ni J m J 

= N n vo\ T 2 n (f- l (tt)) + o(N n ) . 



Here we attempted to apply Lemma 2.5 below as if (/# — z) = for some nice function 
g. As (1.4) shows that is impossible in general. The random perturbation, and taking of 
expected values, make this argument rigorous. In §[4] we show how the first integral split 
to integrals over small (that is of size ~ N~~ 1 ^ 2+e ) subintervals of dil can be replaced by 
integrals of invertible operators. That is done using the singular value decomposition (see 
[18] §3.6] for a simple related example) and facts about random matrices proved in §[3] 
Based on the material reviewed in §[2] we further reduce the analysis to that of traces of 
an inverse of an operator which is a quantization of a slightly exotic function on the torus. 
Here "slightly exotic" refers to the behaviour of derivatives as iV — > oo. An application of 
a semiclassical calculus gives the desired trace and concludes the proof. 

Except for some facts about the standard semiclassical calculus of pseudodifferential 



operator recalled in §2.1[ the paper is meant to be self-contained. One of the advantages of 
Toeplitz quantization is the ease with which traces and determinants can be taken, without 
worries associated with infinite dimensional spaces. 
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random matrices, and Stephane Nonnenmacher for comments on early versions of the paper. 
The authors gratefully acknowledge the partial support by an MU research leave, and NSF 
grants DMS 0500267, DMS 0654436. The first author thanks the Mathematics of U.C. 
Berkeley for its hospitality in spring 2009. 
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2. Quantization of tori 

The Toeplitz quantization of tori, or of more general classes of compact symplectic man- 
ifolds, has a long tradition and we refer to [B] for references in the case of tori, and to [I] 
for the case of compact symplectic Kahler manifolds. We take a lowbrow approach and 
our presentation which follows [15] is self-contained but assumes the knowledge of standard 
semiclassical calculus in IR n . It is reviewed in £2.1 with detailed references to [9] and [TO] 
provided. To see how this fits in the more general scheme see for instance [3 §4.2]. 



2.1. Review of pseudodifferential calculus in IR n . We first recall from [9l Chapter 7] 
(see also [TQl Chapter 3]) the quantization of functions a £ S p (T*M. n ), 

S p (T*R n ) = {a £ C°°(T*M n ) : Va, f3 £ N n , \d^a{x,^)\ < C a ph^ al+mp } , 

< p < - . 
~ H 2 

To any a £ S(T*R n ) we associate its /i-Weyl quantization, that is the operator a w (x, hD) 
acting as follows on ip £ iS(IR n ): 

(2.1) [a"{x,hD)il>Kx) = J J a{^-^)ei^^{y)dyd^ 

This operator is easily seen to have the following mapping properties 

a w (x,hD) : S(R n ) — > <S(R n ) , a w (x,hD) : S'(R n ) — ► S'(R n ) , 

see for instance [TUl §3.1] for basic properties of the Schwartz space S and [TUl §4.3.2] for 
the mapping properties. It can then be shown [9] Lemma 7.8] that a i— > a w (x, hD) can be 
extended to any a £ S p , and that the resulting operator has the same mapping properties. 
Furthermore, a w (x,hD) is a bounded operator on L 2 (IR n ). The condition p < 1/2 is 
crucial for the asymptotic expansion in the the composition formula for pseudodifferential 
operators. If a, b £ S p then 



a w (x, hD) o b w (x, hD) = c w (x, hD), c = a# h b £ S p , 

k 

I z=w=(x,£) j 



°° 1 / h \ k 

(2.2) c{x,0^Y,k\V2^ Dz ' Dw) ) °^ & Hr, 

k=0 ' ^ ' 



where a(z,w) = a(zi, z 2 ,w 1 ,w 2 ) = (z2,w 1 ) - (z 1 ,w 2 ) ■ 



1 fih 



We note that 



so that the expansion in ( 2.2[ ) makes sense asymptotically. 



a(D z ,D w )) a(z)b(w) = 0(h k ^) 
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It is important to recall the standard way in which the quantization of S p (T*M. n ) reduces 
to the quantization of 

S(T*R n ) = S {T*R n ) , 

with a new semiclassical parameter, h = h 1 ^ 2p . Define (x,£) = (h/h)z (x, £), and a unitary 
operator on L 2 (R n ): 

~* n ~ 1 

u h,h u @) = (h/h)*u((h/h)*x) . 

Then 

a(x,hD x ) = U~ia(x,hD £ )U h ~ h , a(z,£) = a((h/h)*(x,£)) . 



We have 



a G S p (T*R n ) -<=>• a G S(T*R n ) 



One simple application of this rescaling is a version of the semiclassical Beals Lemma [9j 
Chapter 8] (see also PH §8.6]): 



(2.3) 



A = a w (x, hD), ae S p (T*R n ) ad^ o • • • ad^A = L 2^ L 2(h N ^) 

for any sequence {£j} =1 of linear functions on T*R n 



The composition formula (2.2) holds also for operators in more general symbol classes. 
For reasons which should become clear below, we will discuss it only for the /i-quantization 
with p = 0. First we need to recall the definition of an order function: m = m(x, £) is an 
order function if there exist C and M such that for all and (£',£')> we nave 

m(x,i) <Cm(x', 0(1 + d R 2 n ((x, £),(£', i'))) M . 

We then say that d G S(rh) if for all a, \d?rd(x,^)\ < C a m(x, £). If rhi and m 2 are two order 

functions and a G S'(mi), 6 G S{fri2) then a(x,hD) o b(x,hD) = c(x,hD), c G ^(mi^), 
and the asymptotic expansion (2.2) is valid in ^(mi^). 

This has a standard application which will be crucial in £|5] 

5eS(m), V(i,0, |a(x,f)| > m(x,f) 

(2.4) „ _ _ 

3h W0 < h < h , a w (x, hD) — b w (x, hD) , b G S(l/m) , 

see for instance [TUl §4. 5, §8. 6]. 

The reason that we presented the order functions on the /i-side is motivated by the fact 
that we need the rescaling of these order functions on the /i-side: we say that m = m(x, £) 
is an /z p -order function if there exist C and M such that for all (x,£) and (x',£'), we have 

(2.5) m(x, < Cm(x', 0(1 + 4- (h~ p (x, £), h~ p (x', £'))) M , 

which means that m(x, f) ^ m (h p x, h p is a standard order function defined above. The 
symbol class is defined analogously, a G S p (m) if d a a = 0(/i~' a ' p m). By the rescaling 
argument the ellipticity statement (2.4) is still applicable if p < 1/2. 
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The following /i p -order function coming from [13, §4] will be essential to our arguments 
here, and in ^5] (Lemma 2.6): 



Lemma 2.1. For a G S(T*R n ) 



is an h 2p -order function in the sense of definition (2.5). In addition, for ip G C^°(1R; [0, 1]) 
equal to 1 on [—1,1], 

(2.6) ( \a(x, Ol 2 + ( ^£P\ ) " e S p (m") . 



h 2 P 



Proof. This follows from the arguments in [13] §4] but for the reader's convenience we 
present an adapted version. We will use the notation (x,£) introduced above, with h = 

h 1 ~ 2p . Let us put F(x,£) = f \a(x,^)\ 2 , so that 771(2;, £) = h 2p m(x,C,) , where 



m\x 



A def 



= h- 2p F(x,Z) + l>l. 



To prove (2.5) we need 
(2.7) 



m{w) < Cm{w'){l + d U 2 n (w,w')) M , 
For |/3| = 1, d^F = 0(h p ^F) and hence 

= -j^d^F = 0(h- p VF) = 0(Vfh) . 

For = 2, d^F = 0(/i 2p ), and hence <9^m = 0(1). By Taylor's formula, 

fh(w') < m(w) + C\j rh(w)d K 2n (w, w') + Cd^n(w,w') 2 
< C(l + m(w))(l + d R 2 n (w, w')) 2 . 



As m > 1 this proves (2.7) with M = 2, and consequently the first part of the lemma. 



For the second part we observe that ip(\a\ 2 /h 2p ) G S p (l ) and hence h 2p ip{\a\ 2 /h 2p ) G 
S p (m). This means that we already have the + case of (2.6). But, 



\a(x, 1 2 + h 2p i, (\a(x, i)\ 2 /h 2p ) > m(x, 0/2 . 



and the — case follows. 



□ 



We remark that by introducing h as a small, eventually fixed, parameter, we can include 
the case of p = 1/2 - see for instance [TJJJ §3.3]. That type of calculus is used in [T3] . 



The last item in this review is a slightly non-standard functional calculus lemma: 
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Lemma 2.2. Suppose that a G S (T*R n ), < p < 1/2, and that ip G C C °°(M). Then 
ip (a w {x, hD)a w {x, hD)*/h 2p ) = q w {x, hD) , q G S p {T*W n ) , 

(2.8) q = go + h 1 " 2 ^ + O s (h°°) , 9j GS„ g (£,0 = ^(K^O)l7^ P ), 
gi (x,0=^(|a(x,Or/^ P )gi(^,0, gie5 p , ^GC C °°(M), #supp^l. 

Proof. This is a simpler version of [T31 Proposition 4.1] which follows the approach to 
functional calculus of pseudodifferential operators based on the Helffer-Sjostrand formula 
for a function of a selfadjoint operator A: 

(2.9) ip(A) = - f(z- A)- l d-Mz)dC(z) , i> e C c °° W , 



7T 



c 



where ^ G C^°(C) is an almost analytic extension of ip, ip\u= ip an d d 2 ip = 0{\ Im^l 00 ) - see 
[9j Chapter 8] and references given there. The reduction to the case given in [HI Theorem 8.7] 
proceeds as follows: the operator a w (x, hD)a w (x, hD)* /h 2p = b w (x, hD), where b G S p (mi), 



where mi is an /i p -order function given by h~ 2p m, where m is given in Lemma 2.1 By the 
rescaling argument above, which gives a reduction to the case of the calculus with h = h l ~ 2p , 
we can apply [91 Theorem 8.7] which gives that ip(a w (x, hD)a w (x, hD)*/h 2p ) = g w (x, hD), 
where g G S p (m^ 1 ) C S p (l). The symbolic expansion presented in P, Chapter 7] complete 
the proof. □ 



2.2. Quantum space associated to T 2n . To define this finite dimensional space we fix 
our notation for the Fourier transform on <S'(M n ): 



= (^72 / < x ) e ~" hM dx . H = ^h 1 , 



and as usual in quantum mechanics, J 7 h u (0 is t ne "momentum representation" of the state 
u. To find the space of states we consider distributions u G <S'(lR n ) which are periodic in 
both position and momentum: 

(2.10) u(x + £) = u(x), F h u{i- + 1) = F h u(£) , 

see [T5| §4.1] and references given there for more general spaces with different Bloch angles. 
Let us denote by TL^ the space of distributions satisfying (2.10). The following lemma is 
easy to prove. 

Lemma 2.3. TC% ^ {0} if and only if h = (2-kN)^ 1 for some positive integer N, in which 
case dim 7^ = N n and 

(2.11) HI = span I -J= 5 ( x ~ 1 ~ H N ) '■ 3 e (Z/N%T 
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For h = (2ttN) 1 , the Fourier transform Th maps Ti 1 ^ to itself. In the above basis, it is 
represented by the discrete Fourier transform 

e -2i-n(j,j')/N 



(2.12) 



N)j,j' 



jyn/2 



j,j'e(Z/NZ) 



The Hilbert space structure on Ti h will be determined (up to a constant) once we define the 
quantization procedure. That will be done by demanding that real functions are quantized 
into self-adjoint operators. 



2.3. Quantization of C°°(T 2n ). The definition (2J) immediately shows that for f e C°° 
satisfying 

V£,m,GZ", f{x + £,ti + m) = f(x,t) => f w {x,hD) : H\ — ► H 1 ^ , 

). Identifying a function / G C°°(T 2n ) with a periodic 



where we consider Ji\ C S' 
function on IR 2n , we define 



f N = f w (x,hD)\ Hl , h 
and we remark that In = Id-n n ■ 



1 



2irN 



C°°(T in ) 3 f 



f N G C(Hl,Hl) 



The composition formula from £2.1 applies since a,b G C°°(T 2n ) can be identified with 
periodic functions on IR 2n ~ T*R. n and 

(2.13) a N ob N = c N , c = a# h b, h = —— : , 

Z7TiV 



where a^^fe is as in (2.2). This means that we simply use the standard pseudodifferential 
calculus but act on a very special finite dimensional space. 

The Hilbert space structure on TC% is determined by the following simple result [T5| 
Lemma 4.3] which we recall below. 

Lemma 2.4. There exists a unique (up to a multiplicative constant) Hilbert structure on 
for which all fx '■ Ti-h ~ > with f & C°°(T 2n ;M.) are self- adjoint. One can choose the 
constant so that the basis in (2.11) is orthonormal. This implies that the Fourier transform 
on TV^ (represented by the unitary matrix {2.12)) is unitary. 



Proof. Let (•, «)o be the inner product for which the basis in (2.11) is orthonormal, and 
put 



Q. 



clef 



l__J26(x-£-j/N) : je(Z/NZy\ 



We write the operator f w (x, hD) on Ti 7 ^ explicitly in that basis using the Fourier expansion 
of its symbol: 

f{x, 0= Yl ^ e 2-«^> + <™,0) . 
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For that let L irn (x,^) = (£,x) + (m, £), so that 

f w (x,hD)= J2 f(e,m)exp{2mL^(x,hD)). 

We also check that 

exp (2iriLl m (x, hD)) Q 3 = exp(7u(2(j, £) - (m, l))/N) Qj- m , 
(note that j G (Z/NZ) n and j — m is meant mod iV) and consequently, 

m£Z n /(NZ) n 

F mj = Yl ^ ~ m - rN)(-l)M exp(™ (j + m,£)/N) . 



Since 



Fjm= f(-lj-™ + rN)(-l)M exp(-7ri(j + m,£)/N) 
e,rez n 

exp(7u (j +m,£)/N) , 

we see that for real /,/ = /, Fj m = F m j. This means that f w (x,hD) is self-adjoint for 
the inner product (•, •) . We also see that the map / i— > (i ? J ' m ).j,me(z/-ft/'Z) n is onto, from 
C°°(T 2n ; M) to the space of Hermitian matrices. 

Any other metric on TC% could be written as (u,v) = (Bu,v)o = (u,Bv)o. If (f^UjV) = 
(u, fN v ) for all real /'s, then Bf^ = f^B for all such /'s, and hence for all Hermitian 
matrices. That shows that B = eld, as claimed. □ 



We normalize the inner product so that the basis specified in (2.11 ) is orthonormal. From 
now on we use this basis to identify 

HI - £ 2 {Z n N ) ~ C Nn , Z n N = {Z/NZ) n . 



The calculation of the matrix coefficients in the proof of Lemma 2A immediately gives 
the following 



12 T.J. CHRISTIANSEN AND M. ZWORSKI 

Lemma 2.5. Suppose f G C°°(T 2 ™). Then 

tr f N = N n (-l) mm) fm,Nm) 

£,m£Z n 

(2.14) =N n [ f(w)dC(w)+r N , 

\r N \<C kn N- k+n V / \d a f(w)\dC(w), 

|a|<max(fc,2n+l) 

for any k. Here C(w) is the Lebesgue measure on T 2n normalized so that £(T 2n ) = 1. 

It is well known that for / G C°°(T 2n ), independent of N, Jn is uniformly bounded on 
£ 2 (Z^r) - see [6]. We will recall a slight generalization of that for functions which are allowed 



to depend on N in a S p -w&y described in £2.1 



2.4. S p classes for the torus. The S p classes for the quantization of the torus have 
already been considered in [16] and we refer to that paper for more detailed results such 
as the sharp Garding inequality. Here we continue with a self-contained presentation. 

We first define a class of order functions: a function of w G T 2n and a > is an a-order 
function if there exist C and M (independent of a) such that 

(2.15) V^'GT 2 ", m( ; W ; a \ < C(l + d T n (w/a, w'/a)) M , T™ = f (R 2n /(Z/a) 2n ) , 

with the distance induced from the Euclidean distance: d^2n/ r (w, w') = inf 7g r \w — w' 
With this definition we have 

(2.16) S(m,a) = {aeC°°(T 2n ), d p a(w) = 0( a -^m(w,a))} . 

If 

(2.17) N- p /C <a <CN~ P , 0<p<^, 



the quantization procedure described in {2.3 applies to S(m,a): we now quantize func- 



tions / which are periodic and belong to S p with h = 1/(2ttN). Similarly, we have the 



composition formula (2.13) with the asymptotic expansion in (2.2) valued in S(mirri2, a 



Lemma 2.6 translates into this setting and will be used in 



Lemma 2.6. For f G C°°(T 2n ) 



m(w, a) = \f(w)\ 2 + a 2 , 
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as an a -order function in the sense of definition (2.15). In addition, for ip G C^°(IR; [0, 1]) 
equal to 1 on [—1,1], 



(2.18) 



l/ff + ^l^feS^a). 



For 5(1, a) we also have uniform £ 2 -boundedness, which we present in the simplest form: 



Proposition 2.7. Suppose f G 5(1, a) with a satisfying (2.17). Then 
(2.19) \\f N \\p^P < sup |/| + o(l) , iV^oo. 

f2n 



Proof. Lemma 2.5 gives 



II f II 2 def 

II/amIhs ~~ 



tif* N f N = N n [ f# h fd£ + 0(N~ k+n ) V / \df>(f# h f)\dC, k^n. 
Jt*» \0\< k Jr2n 



Since fj^hf G 5(1, a) (that is, using (2.17), /#/i/ lies in S p when considered as a periodic 
function on IR 2n ), we see that 

H/jvllls = 0(jV n ) + C(A^- fc(1 - p)+n ) = 0{N n ) . 

Hence, 

(2.20) ||/iv||^<||/iv||HS<CiVf. 

We now use Hormander's trick for deriving L 2 -boundedness from the semiclassical calculus. 
Let M > sup T2 „ |/# & /| and let a N = M - f* N f N , a = M- f# h f G 5(1, a), a > 1/C > 0. 
Then by Q 

b° N b% - a N = r% , r°GAT 2 ^ 1 5(l,a), 6° == \/a G 5(1, 



We now proceed by induction to construct real V G N^ 2p ~^S(l, a), < j < J, so that 

(2.21) (i^) 2 - = ^ , l^ d = f X> w , r J GiV^ +1 )^- 1 )5(l,a). 

i=o 

Suppose that we already have it for J (the first inductive step being J = 0) and we want 
to find b J+1 G A^( J+1 )( 2p - 1 )5(l,a) so that : 

(B J N + ^ +1 ) 2 - a N = r J N + B J N b J N +1 + b J N +l B J N 

= 4 + b° N b J N +1 + b J N +1 b° N + R J N b J N +1 + b J N +1 R J N , 

where R J = B J — b G A r2p_1 5(l, a). We now simply put 

b J+i = r J/(2b°) G iV (J+1)(2 ^ 1) 5(l,a) , 



which is real since the left hand side of (2.21) is self-adjoint. The inductive step follows 
again from the composition property. 
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Returning to the boundedness on C 2 we now have 

M IM| 2 - ||/jv«|| 2 = (a N u,u) = (B J N u,B J N u) - (r J N u,u) 
> -\\rjr\\p-+p \\u\\ 2 > -|l r )vllHs|M| 2 

>_ C r JV S+(J+l)(2p-l)|| w ||2 



where for the last inequality we used (2.20). Hence by taking J large enough, ||/jv||^ 2 -»f 2 < 
M 1//2 + o(l), and since M can be taken as close to sup |/| as we like, this gives (2.19) □ 

One of the consequences of the boundedness on I 2 is the justification of the basic principle 
of semiclassical quantization: 

Poisson brackets, {•, •} < — ► Commutators, (i/h)[» N ,» N ], h = l/(2irN). 

More precisely, {/, g} = ^2] =1 {d^fd Xj g-d x .fd^g), (with o = J2] =1 d^Adxj the symplectic 
form on T 2n ), and 

(2.22) 27riN[f N ,g N ] = ({f,g}) N + Op^(N~ 2+A P) . 

The functional calculus lemma presented in the W 1 setting translates to the case of the 
torus: 

Lemma 2.8. Suppose f G C°°(T 2n ) and a = h p , < p < 1 /2. Then, for G C C °°(M) ; 

(2 - 23) q = q° + ^- 2 V + Os(h°°) , q j eS(l,a), q°(w) = ^(|/H| 2 /a 2 ) , 

q\w) = m(w)\ 2 /a 2 )q\w), ?eS(l,a), $ G C C °°(R) , #su PP ^ 1 • 

Proof. We need to check that for a function <p G C£°(R), and g G C°°(T 2n ; R), the action 
of (p(qn) 011 — ^ 2 (^jv) defined using functional calculus of self-adjoint matrices is the 
same as the action of (p(g w (x, hD)) on Ti 1 ^ C iS'(IR n ). In view of the Helffer-Sjostrand 
formula that follows from verifying that the action of the resolvent (z — g^v) -1 ; Imz 7^ 0, 
on Ti 1 ^ is the same as the action of (z — g w (x,hD))~ 1 , Imz 7^ 0, on Ti 7 ^ as a subset of 
5(lR n ). But we know from (2.3) that for Imz 7^ 0, (z — g w (x, hD))~ l = F(z,x, hD), where 
F(z) G 5(1) (non-uniformly as Imz — > but with seminorms polynomially bounded). 
This means that the L 2 inverse is a restriction of an inverse defined on S'(M. n ). Hence 
(z — gw) 1 = [F{z)]n and the actions are the same. This argument is not asymptotic in N 
and applies to ip — ip{»/a 2 ) and g = fj^hf- D 



Proposition 2.9. Suppose that (1.2) holds with z — 0. Then for any ip G C£°(R), 
(2.24) rank^ ( < CN n a K , N~ p < a < 1 , p < 1 



a J 2 
wfflt/i i/ie constant depending only on the support ofip. 
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We note that by proceeding either as in the proof of [131 Proposition 4.4] or as in the 
proof of [121 Proposition 5.10] we can show that the result is valid for p = 1/2 but we do 
not need that in this paper. 

Proof. Suppose ipi G C^((-R 2 + 1,R 2 - 1), [0, 1]), R » 1 is equal to 1 on the support of 
ip. Then, using the functional calculus of self-adjoint matrices and Lemmas 2.5 2.8 , and 
(B we get, with $ G C c °° ((-#*, R 2 ), [0, 1]), # S u PP ^= 1, 

rank ^M^ < tr ^i (^^j < N n J ^(\f\ 2 /a 2 )d£ + 0(N~°°) 

< N n £({w : \ f(w)\ < Ra}) + 0(N-°°) < CN n a K , 
proving the lemma. □ 



3. Some facts about random matrices 

Random matrix theory is a very active field and we refer to Mehta's classic book [H] 
for general background, and to [H] for some recent works and applications. All the facts 
we need in this paper are elementary but they do not seem directly present in the main- 
stream literature. Consequently the presentation is almost self-contained and, reflecting 
the authors's own position, does not assume any knowledge of the subject. 

We consider the ensemble of complex Gaussian matrices with independent entries dis- 
tributed in C according to the standard normal distribution. That means that there exists 
a probability space, (Q, E,/i), Ea cr-algebra of subsets of Q and //:£—> [0, oo), a measure, 
with yu(fi) = 1, and a map Q 3 uj i— > Aj^(u), An{uj) = {aij{uj))i<ij<di such that uj \— > a,ij(cu) 
are independent random variables with standard normal distribution. The push forward 
measures on C, (ajj)*/i, are given by exp(— \z\ 2 )d£(z)/n, where £ is the Lebesgue measure 
(standard normal distribution), and 

[(aii,aM)].M = \e-\ z \ 2 -\ w \ 2 d£(z)d£(w), + (k,£), 

(ciij, au) : Q — > C z x <C W , which is the statement that a%j and are independent. 

A more useful global description of the random variable Ad{oj) is given as follows: let 
= (an, ...,a id y G C d , and set A = (ai, a^). Denote 

d 

d£(ai) = dHeandIma i i...dIieai C idlmaid, and d£{A) = Y\_d£(a,i). 

i=i 

Then, as a measure on C d \ the space of d x d matrices, 

(3.1) A*n = 7c- d2 exp (-||A||^s) d£{A) , ||A||| S = tr A* A , 
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3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 



Figure 3. A numerical example suggesting that Proposition 3A_ is optimal: 
the left hand side is computed numerically for A = diag ([0,0, rand (1,8)]) 
(a 10 x 10 diagonal matrix of rank 8) where rand command produces uniform 
distribution on [0, 1]. It is plotted as a function of log(l/<5). The upper bound 
in Proposition 3.1 (with C — 1) is also plotted for comparison. 

where HS stands for Hilbert-Schmidt. Note that each entry of A is a complex iV(0, 1) 
random variable. 

We recall that any matrix A can be written using its singular value decomposition, 
(3.2) A = USV* , 

where UU* = U*U = Id, VV* = V*V = Id, that is U and V are unitary, and S is a 
diagonal matrix with non-negative entries. If the entries of S are distinct and we order 
them, the decomposition is unique. 

Proposition 3.1. Let A be a constant d x d matrix, and let Q be a d x d random matrix, 
with the entries qij independent complex N(0, 1) random variables. Then there exists a 
constant C independent of d and A, such that 

J |E(tr(L4 + 5Q)~ 1 A) \dt<Cti log ( 2 + ^f) ) ' 

where \A\ = ^/AA*. 

In the proof we will need the following 
Lemma 3.2. The function g(s) = f j c (l/\s + q\)e~^ 2 dC(q) is smooth for s6C, and 

g(s) = j^j- + O ^j~|2^ ' as \ s \ ~ > °°- 

Proof. The asymptotic expansion follows from the integrability of |g|, a change of variables, 
q i— > 1 + q/s, and the method of stationary phase. □ 
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Proof of Proposition 3.1 Using the singular value decomposition for A, we may write A = 
USV*, with U, V unitary and S a diagonal matrix with non-negative entries a\, cr^ on 
the diagonal. We note that 

tr((L4 + 5QY 1 A) = tv((tUSV* + 5QY l USV*) = tr((tS + 5U*QV)- l S). 

Since U*QV is a random matrix with the same probability distribution function as Q, we 
have 

E(tr((M + 5Q)- l A)) = E(tr((tS + SQ^S)). 
Thus we may assume that A is diagonal, with non-negative entries <Ti, cr^. We have 

where here and below Mij is the (i, j) minor of the matrix tA + 5Q. 
To compute E(Mjj<7i/det(ty4 + 5Q)), we write 

det(M + 5Q) = (ten + 5 qii )M u + ^(-1)^%M, 



and define 
(3.3) 



clef 



qeC d2 : {{tcTi + Sq^Mal > 



Let If be the characteristic function of a set F. Then 



(3.4) 



E 



(- 



Mao; 



E 



Mad 



\det(tA + 5Q) 
since the boundary of has measure 
Now, 



det(L4 + <5Q) 



11 



\det(tA + 5Q) s * 



E 



M,;,;<7,; 



det(M + 5Q) 



h 



E 



E 



MnOi ( E-^(-i) i+ % 



oo 

E 



fc=0 



(to-j + SqujMu 



■l-This follows from the fact that the pushforward of the probability measure by Q (the probability 
density) is absolutely continuous with respect to the Lebesgue measure on C™ and the set 

d 

{QeC" 2 : Q = («tf)i< W <n, |(tffi + a<te)M«| a = ^(-ly+^^-Myl 2 }, 
has Lebesgue measure 0. 
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We recall that the set E& is chosen so that the infinite sum converges. 

The set E& is invariant under the mapping 

(3.5) q a , g M _i, q i>i+ i, g i)d i-> e^ga, e^g M _i, e^g^+i, e v g ijd 

for any real number yj. Since My's are independent of g^, Ej^i( — * s homo- 
geneous of degree 1 under this same mapping and (to~i + 5qu)Mu is independent of g^- for 
j 7^ z, we find that 



E 



det(L4 + 5Q) 



E 



(t(Ti + 



-1 



We do a similar computation for the second term of (3.4): 

(- 



Vdet(tA + 5Q) E ' 



E 



E 



Mad 



MaOi 



(toi + 5qu) 



-i 



E^(-i) i+ ^%^- 



E-^(-i) i+ ^%^, 



Is 



using, as before, the invariance properties of and the homogeneity of 
Thus we have 

d 

(3.6) E(tr(tA + 5Q) _1 A) = ^E 

Now, 



0. 



i=i 



Or 



(t<7i + 8qi 



E 



(t<Ti + 



-I 



"^•11 







j dt 









Ti/S 



E 



toi + 5gjj 
1 

|s + qu\ 



dt 



ds 



E 



71 



Vl t0 "iA + fel 

cr i; /5 

g(s)ds 



dt 



where g is the function defined in Lemma 3.2 Using this, (3.6), and the results of Lemma 
3.2 proves the proposition. □ 



Lemma 3.3. Let F, G be d x d matrices, with F invertible, and let (3 = \\F 1 ||. Then 
E(tr ((F + SQ^G)} = tr (f^C?) (l + O^ 1 / 4 ^ 2 )) + O ^G^e' 1 '^ 5 ^ . 
T7z,e implicit constant in the error term is independent of F and G. 
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Proof. We first note that if we replace F by its singular value decomposition, F = USV*, 
then 

E(tr ((F + 5Q)- l G)^j = E(tr ((S + 8Q)- 1 {U* GV)f) 

and 

tr (F^G) =tr (S^ITGV). 
Thus we may assume that F is a diagonal matrix. 



Our proof then resembles the proof of Proposition 3.1 Let x £ be the charac- 

teristic function of (—00, 1/2], and, if A = (ay), let ||A|| sup = sup^- |ay|. We write 



(3.7) E(tr((F + 5Q)- 1 G))=E(tr((F + 5Qy 1 G)x(d\\Q\\ snp 5f3) 

+ E(tr ((F + ^g)- 1 ^) (1 - xMIQIIbup^)))- 
For the first term, 

00 

E(tr((F + <5g)- 1 G)x(rf||g||sup^)) = E(tr (^^(-^O^^xCdligiU^)))- 



Using the fact that the cut-off x(d||Q||sup<fy3) is invariant under rotations of the gy and that 
the qij are complex and independent, we find 



E(ta ((F + 5Qy 1 G)x(d\\Q\U P 5P)) = tr (F^GMQ : ||Q|| sup < l/28(3d) 
(3.8) = tr (F- 1 G)(1 + 0(dV 1/4(w2 )) . 



Now we consider the remaining term of (3.7). In a way similar to the proof of Proposition 



3.1 we denote the diagonal entries of F by fa = and by My the minor of F + 8Q. 



If G — (gij), we have 

;(tr {(F + 5g)- x G)(l - X (d\\Q\\ sup 6(3))) = ^( ^aeUF + SO) ^ ~ x(*IU*0))- 



El 



h3 



Just as in the proof of Proposition |3.1[ to compute 

Kd^) (1 -* (d|IOIUW) ) 

we write 

det(F + SQ) = fa + 5 qii )M u + ^(-l)^yMy 



and define Ejj as in (3.3). Proceeding almost exactly as in the proof of Proposition 3.1 



using that both £ i4 and the support of (1 — x(d\\Q\\ sup 5 (3)) are invariant under the mapping 
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(3.5 ), we get that 

Mugu 



E 



But 



det(F + 6Q) 



E 



[l-x(d\\Q\\ sup 6f3)) =E 



fa + Squ) 



9i 



{(Ti + Sqi 



-lx u (l - X (d\\Q\\ sup 5(3)) 



Wl-xWIU^)) ■ 



I CI 



o 



To compute 



when % ^ j, we write 

(3.9) det(F + 5Q) = 5 qji M ri (-l) i+j + (a, + Squ) Ma + £ (-l) fc+i 5 fc M, 
and define 



fct 



def 



g G C* 2 : |%M^| > 



-'j i 



i^k^j 



ki 



Following the proof of Proposition 3.1 but treating the term bq^Mji as the distinguished one 
in the expansion of the determinant (3.9) and using the invariance of £jj under rotations 
of qji, we find that 



( {-!)*+* Mji9ji 
Vdet(F + 5Q) 1 



X (d\\Q\U P SP)) 



E 



{a { + Sq u )M u + Ek&A-V^QkiMki 
Since on the support of % c 

1 



l^(l- x (d\\Q\\ mp Si3))\. 



Mji\ < 



a i + Sq ii )M ii + J](-l) fc+ Mg fci M; 

k^i,j 



ki i 



we find 



E 







□ 



Our proof of Proposition 4.1 in the next section will use Proposition 3.5 To prove this 
proposition we will need several preliminary results. 
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The first lemma below follows from well-known facts about eigenvalues of complex Gauss- 
ian ensemble. We give a direct proof suggested to us by Mark Rudelson: 

Lemma 3.4. Let A = (a\, ...,a d ), with ai G C d . Then, with the notation of (3.1), 

/ \ det Al' 1 d£(A) <oo. 

Proof. We begin by introducing some more notation. For p < d, p G N, v G C d denote by 
V p v projection onto the subspace spanned (over the complex numbers) by ai,...,a p . This 
of course depends on a 1; ...,a p , but we omit this in our notation for simplicity. 

Using the Graham-Schmidt process, we can, if A is invertible (as it is off a set of measure 
0), write the matrix A = UR, with U a unitary matrix and R being upper triangular. The 
diagonal entries of R are then given by ||ai|| and ||(1 — V p -i)a p \\, p = 2, ...,d. Thus 

| det A | = ||ai||||(l-Pi)a 2 ||||(l-P 2 )a 3 || •■■ ||(1 - P d -i)a d \\. 

Note that 

IKIIIKl -Pi)a 2 ||||(l -V 2 )a 3 \\ ■ ■ • ||(1 -V d -2)a d -x\\ 
is independent of a d , that is, independent of aid, d2d, ■ ■■,a> dd - Therefore 

/ \&etA\- l dC{A) 

J\\A\\hs<1 

1 dC(a d )d£(a d -.i)...d£(ai) 



< 



\\A\\ H s<i I^iIIIK 1 _ ^i) a 2|| • ■ • ||(1 - V d -2)a d -i\\ 11(1 - V d -i)a d \\ 

f dC{ad) dC{a d -i)...dC{ai) 

\\ax\\<i J\\a d \\<i 11(1 -Vd-i)a d \\ ||ai||||(l -Vi)a 2 \\ ■ ■ ■ ||(1 - 7V 2 )a d _i| 



The value of L ,, K1 1/||(1 — V d ~i)(id\\dC(ad) depends only on d and the rank of the space 

spanned by ai,...,a d -i- We find 1/||(1 — V d -i)a d \\ is locally integrable over M. 2d ~ C d , 
because G C d and the space spanned by ai,...,a d -i has complex dimension at most 
d — 1. Therefore 

(3.10) f 1 dC(a d ) <C<oo. 

J\\a d \\<i IK 1 ~ r d -i)a d \\ 

Here the constant C can be chosen independent of ai, a d -i, as the maximum of the 
integral in (3.10) occurs when ax, a^_i span ad-l dimensional vector space. The proof 
follows by iterating the above argument. □ 

Proposition 3.5. Let A(s,t) be a dxd matrix depending smoothly on (s,t) G U C C 2 . Let 
Q denote a d x d random matrix, with each entry an independent complex N(0, 1) random 
variable. Then for 5 > 0, (s,t) G U, E(tr((A(s, t) + 5Q)~ 1 d t A) is smooth on U , and 

d s E(tr{{A{s,t) + 5Q)- 1 d t A)) = d t E(tx{{A{s,t) + SQY^A)). 
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This proposition has the following corollary. 
Corollary 3.6. Let M , B, be d x d matrices independent of s and t. Then 

J E^tr ((sB + M + 8QY l B)}ds = E ^ tr ((B + tM + 5Q)~ l M) ^ dt 

-J E^tr ({tM + 5Q)- 1 M)^dt + J E^tr ((sB + 8Q)~ l B)^jds. 

Proof. Using the previous proposition, this follows from the Fundamental Theorem of Cal- 
culus: 

J E^tr ((sB + M + 5Q)- 1 B)^)ds - J E^tr ({sB + SQ^B^ds 

= [ d t [ Eti ((sB + tM + 5Q)- 1 B) dsdt 
Jo Jo 
1 ,i 

d s Etr ((sB + tM + 5Q)~ 1 M)dtds 
Jo 

E^tr ((B + tM + SQ^M^dt -J E^tr ((tM + SQ^Mfjdt. 



□ 

Proposition 3.5| follows from the subsequent two lemmas. 



Lemma 3.7. Let A(s,t), B(s,t) be d x d matrices depending smoothly on (s,t) G U C C 



2 



With Q a random matrix as in Proposition 3.5 and S > 0, 



E(tr ((A(s,t) + 5Q)- 1 B(s,t))) e C°°{U). 

Proof. We prove the lemma by writing the expected value as an integral: 
E(tr((A + 5g)- 1 fi)) = J tr((A + 5Q)- 1 B)e-^HsdC(Q) 



J ti((5Q)- 1 B)e^ Q - 1 s A ^sdC(Q). 



Now, for a d x d matrix B, | tr^Q)" 1 !?)! < C\ detQ| -1 ||S||||Q|| << - 1 /^ where the constant 
C depends on d. Moreover, 

j+k 

j'<j,k'<k 

I 4112 



\did^-^H<c jlk>d ( £ ii9r^Aii)(M) J+fc e - 

.V.' n<-t \ / 



Since, using Lemma 



3.4 



/ 1 det + \\Q\\) m e-\\ Q -s A W»sdC(Q) < oo, for any finite m, 



the smoothness of A and B proves the lemma. □ 
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If M is an invertible matrix depending smoothly on s and t, then 



(3.11) 



/ w 1 o , ^ <9t det M 



and d s ti(M~ l M t ) = d t ti(M~ 1 M s 



detM 

The lemma below shows that something similar is true when taking expected values, even 
though the matrices under consideration are not invertible for some values of the random 
variable. 

Lemma 3.8. Let A(s,t) be a d x d matrix depending smoothly on (s,t) G U C C 2 ; and Q 
a random matrix as in Proposition 3.5\ Then for 5 > 

d s E (tr((A + 5Q)- l d t A)) = d t E (tr((A + SQy'dsA)) . 

Proof. Let Xe £ C°°(M) satisfy Xe( x ) = 1 f° r l x l < e /2 and Xe( x ) — for |x| > e. Then 

d s E (tr((A + 8Q)- l d t A)) = d s E ( X e(det(A + SQ)) tx((A + SQY 1 d t A)) 
(3.12) + 8 S E ((1 - x,(det(A + 8Q))) tr [(A + SQ)~ 1 d t A)) . 

Now 

d s E ((1 - X e(det(A + SQ))) tr ((A + SQY^A)) 

(1 - X e(det(A + 6Q)))d 8 tr [(A + SQ)~ l d t A)e- m »sdC{Q) 



x' e (det(A + SQ)) (d a det{A + SQ)) tr ((A + SQ)-%A)e-^BsdjC(Q) 

where we can freely interchange differentiation and integration since the integrand is smooth 
and it and its derivatives are integrable. But using (3.11), we get 

8 S E ((1 - x,(det(A + SQ))) tr ((A + SQ)~ 1 d t A)) 

(1 - X e(det(A + SQ)))d t tr ((A + SQY 1 d s A)e~ m ^dC(Q) 

- J x' e (det(A + 5Q))d t det (A + SQ) tr ((A + SQ)- l d s A)e^ Q ^sdC(Q) 
= 8 t E ((1 - Xe{det{A + SQ))) tr ((A + dQ^^A)) . 

On the other hand, the first term on the right in ( 3.12[ ) satisfies 
lim<9 s E (xe(det(A + SQ))(tr((A + SQ)~ 1 d t A)) 

= \imd s [ Xe(&A(A + 5Q)){tr((A + 6Q)-%A)e-^BsdC(Q) 
40 J 

= limd s [ xMet(SQ))(ti((SQ)- 1 d t A)e^ Q - 1 ^sdC(Q) = 
40 J 

since (tr((<5Q) _1 d t A)e~^~ '5 A "« s and its s derivative are both in L 1 , using Lemma 



3.4 



□ 
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4. Reduction to a deterministic problem 



In this section we will show how to reduce the random problem problem to a deterministic 
one. That will be done using the singular value decomposition of the matrix f^. 

Let A be a square matrix, and let USV* be a singular value decomposition for A. We 
make the following simple observation: for ip G C£°(R, R) equal to 1 on [—1, 1], 

(4.1) (A + aiPiAAya^UV*)- 1 = 0(1/ a) : i 2 — ► f , 

which becomes totally transparent by writing tp(AA* /a 2 )UV* = Uip((S/a) 2 )V*. 

The random problem is reduced to a deterministic one by using an operator of the form 
(0). 

Proposition 4.1. For a smooth curve 7 define 

(4.2) I N (<y) = j Etr(/jv + SQn ~ z)~ l dz 

where Qn is a complex N n x N n matrix, with entries indepent N(0, 1) random variables. 
Let /a? = UnSnV/j be a singular value decomposition of fx, and let ip G C£°(R; [0, 1]) be 
equal to 1 on [—1,1]. // 

(4.3) G 7, |7| < a/4, 5 < a, 
then 

/ Etr(/jv + SQ N (u) - z)~Hz = 

J-y 

(4.4) / EtrC/W + ai,{f N r N /a 2 )U N V* N + 5Q N (u) - z)~ x dz + E l = 
tr(f N + a^(f N f* N /a 2 )U N V^ - z^dz + E 2 , 



®\ _|_ r -a 2 /4(37V"<5) 2 
.5 



w/iere 

(4.5) £ 1 ,£ 2 = o(aTog(^) + V e 

and d = rank ]l supp ^ (/n/n/® 2 ). 



The proof of this proposition will use the following lemma. 
Lemma 4.2. Let fx, Un, Sn, Vn, if), 5, d, and a be as in the statement of Proposition 



4-1 Let x G L°°(R) be the characteristic function for the support ofip. Then, if \z\ < a/4, 
1 E tr ((/jv + sa X (f N r N /a 2 )U N V* - z + 5QY 1 a X (f n T N / a 2 )U N V* N ) ds 



satisfies the bound (4.5). 
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Proof. First suppose that for a m x m matrix A, 

(4.6) A=(f f ) and J- = ( f' f> 

V ^21 ^22 / V -021 -D22 

with An, Bu d x d matrices and An, Bn (m — d) x (m — <i) matrices. Then if A22 is 
invertible, we have the Schur complement formula, 

(4.7) B n = (An + A 12 A 22 l A 2 i 
see [18] for a review of some of its applications in spectral theory. 

We note, using ip(AA* /a 2 )UV* = U^((S/a) 2 )V* and the unitarity of U N , V N , 

(4.8) Etr ((f N + sax(f N f* N /(y 2 )U N V* - z + 5Q N y l axUn F N I a 2 )U N V N ) 

= Etr ((S N + sax{S N S* N /a 2 ) - U* N zV N + 5Q N )~ l ax(S N S* N /a 2 )). 

The main idea of the proof will be to effectively reduce the dimension of the matrices we 
work with, from N n to d. We can assume that Un, Vn, Sn are chosen so that the diagonal 
elements ax, &N n of Sn satisfy <7i < cr 2 ■ •• < a N n. Let J denote projection onto the 
range of x(5jv/a 2 ), which is the same as projection off of the kernel of x(S%/a 2 ). Then 

Id 




J 

and ax{S N /ot 2 ) takes the form 



We also write 



al d 




C 1 / C C* / 2\ tt* \r ( s<y Id + An A12 

S N + sax{S N S N /a ) - U N zV N = I ^ ^ 



and 

\ Q21 Q22 

where An, Qn are <i x d- dimensional matrices, and A22, Q22 are (N n — d) X (iV n — <f)- 
dimensional. Since is diagonal and \z\ < a/4, we have ||A 12 || < a/4, ||A 2 i|| < a/4. 

Using this notation, we have that A 2 2 is invertible, with norm at most 4/3a. Now restrict 
Qn to the set with 

(4.9) S\\Q N - JQ N J\\ sup < aN- n /A. 

Note that poses no restriction on Q n . For such Q N , A 2 2 + 8Q22 is invertible, with norm at 
most 2/a. Restricting to this set of Qn and using (4.7), we find 

tr US N + sax{S N S* N /<y 2 ) - U N zV N )~ l ax{S N S N l 'a 2 )) = tr d (a{sal d + M d + SQn)' 1 ) , 
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where we use the notation tr^ to emphasize we are taking the trace of a d x d matrix, and 
where 

M d = A u - (A 12 + 5Q 12 )(A 22 + 5Q 22 )-\A 2l + 5Q 2l ) 

is a d x d matrix depending on Qi 2 , Q 2 i, and Q 22 , but not on Q n . Since ||Au|| = {^(Sn — 
zU n Vn)J\\ < Ca and \\A\\i 2 < a/4, ||v4|| 2 i < a/4, we have \\M d \\ < Ca, for a new constant 
C independent of N, 6, and Qn satisfying (4.9). 

Next we take the expected value in the Qu variables only: 

E Qll (F(Q N )) = [ F(Q N )e^ Q ^sdC(Q n ). 

Still requiring Q N to satisfy (4.9), which is not a restriction on Q u , and using Corollary 
3.6[ we get 



J E Qll (atr d (M d + saI d + 6Q u y 1 ))ds = 
J EQ^tu^Md + ah + dQ^YHl^dt- E Qll ^tv d (saI d + 8Q ll )- 1 a)y s 
+ ^ E Qll (tr d (tM d + 6Q 11 )- 1 M d )yt. 



Recalling that \\M d \\ < Ca we see from Proposition 3.1 that the second and third terms 
on the right are 0(dlog(a/8)), if a/5 > e. Moreover, 



\M d -JS N J\\ < 



a 



and Sn > 0. Therefore, for < t < 1, al d + tM d is invertible, with the inverse having 
norm at most 3/a. Thus from Lemma 3.3 we see that 



J E Qll (tr d (tM d + al d + 5Q 11 y 1 M d fjdt 



0(d) + O ( ^e"" 2 / 4 ^ 



The implicit constants in both cases are independent of Q — JQJ satisfying (4.9). Thus 
we get 



(4.10) 



E(tr ((S N + sax(S 2 N /a 2 ) + 5Q - zU* N V N ) l ax{S 2 N /a 2 ))t {m ^j Q j\ lmp <^ w} 



ds 



= 0(dlog(a/8)) + O [d i 5- l e- a /4(d35y 
where for a set E, lg is the characteristic function of E. 
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Exactly as in the proof of the Lemma 3.3, we can show that 

e( tr ((S N + sax(S 2 N a 2 ) + 5Q N - zU^V^axiS^a^t^Q^jQj^^/^n^ 

= O (V 4n <r 1 e" Q2/(4Arn<5) 



(4.11) 



Using (4.8), (4.10), and (4.11), we prove the lemma 



□ 



We now use Lemma 4^ in a preliminary step towards proving Proposition 4.1| 
Lemma 4.3. Let /at, Z7 at, Sn, Vn, ip, 5, d, and a be as in the statement of Proposition 



4-1, and set x = % 



supp tp ■ 



Then 



Etr(/W + SQ N - z)- l dz = / Etr(/jv + a X (f n f* N / a 2 )U N V* N + 5Q N - z)~ l dz 

_/V 4ra 



+ o(dlog(|))+0 



-a 2 /A(3N n Sy 



5 



Proof. The proof uses the same type of argument as Corollary |3.6| Using the Fundamental 
Theorem of Calculus, 

/ EtrC/W + 5Q N - zY l dz - [ Etr(/jv + u X {f n f* N / a 2 )U N V* N + SQ N - z)- l dz 
= - / d s I Eix(f N + sa X (fNfN/<x i )UNVZ + 6Q N - Z )- 1 dzds 

JO J-y 

= - [ d z [ E (tr(f N + sax{f N r N /a 2 )U N V; + 5Q N - z)' 1 a X (J n f* N / 'a 2 )) dsdz 

J-y JO 



where we use Proposition 3.5 The right hand side is 



Y> / Eti ((f N + sa X (f N r N /a 2 )U N V^ + 6Q N - z^axifNf^/a 2 )) ds 
± Jo 



where z± are the endpoints of 7. Then using Lemma 4.2 finishes the proof. 



□ 



We are now able to give a straightforward proof of Proposition 4.1 



Proof of Proposition ^. 1\ We begin by noting that, with x — Hsu PP »/> 

WUn + ax(f N r N /a 2 )U N V* - z)' 1 ]] = 0(1/ a) 

and 

\\(f N + a^(fNr N /a 2 )U N V* N - z)- l \\ = 0(l/a) 
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when \z\ < a/4. Moreover, the rank of x(/Af/]v/ a2 ) i s ^ an d the ran k of ^(/Ar/iv/^ 2 ) * s a ^ 
most d, and both operators have norm at most 1. Then 



(tr(/jv + ax(fNf* N /a 2 )U N V^ - z)' 1 - tr(f N + atfj(f N f* N /a 2 )U N V^ - z)' 1 ) dz 

tr ((f N + a X (f N r N /a 2 )U N V* - z)" 1 (x(f N F N /a 2 ) - ^(f N F N /a 2 )) 
x (f N + a^(f N f* N /a 2 )U N V^ - zY 1 ) dz\ 

< [ 9± dz = o(d). 



Thus, applying Lemmas 4.3 and 3.3 proves the Proposition. 



□ 



5. Proof of Theorem 

The proof of Theorem will be deduced from the following local result: 

Proposition 5.1. Under the assumption of the main theorem, let 7 C dQ be a connected 
segment of length 



(5.1) 



a 



1 ,n 1 

a = h p , < p < - 



and let In(i) be as defined by (4-2). Then for exp (— h e ) < 5 < h Po , we have 



(5.2) I N (j) = N n (f(w) - z)- l dC{w)dz + 0{\-f\h- n+f) ^- 2t ) + 0(|7|/r n+1 - 2 O 

Jy Jf2n 



where we note that (1.2) with k > 1 implies that (f(w) — z) 1 G L l {T 2n ) so that the first 



term on the right hand side makes sense. 



Assuming the proposition we easily give the 



Proof of Theorem. We divide dQ into J = C /a disjoint segments jj, < a/C. Propo- 
sition 5.1| implies that 



N 1 



E ( tr Ln fN + 6QN ~ Z)1 ) = J^ In ^ = 

[ [ (f(w) - z)- l dC(w)dz + 0{h- n+p ( K -V- 2e ) + 0{hT n+l ~ 2p ) 
Jon JJ 2n 



We now choose p = 1/(k + 1), to optimize the error, that is to arrange, p(n — 1) = 1 — 2p. 
That means that the error is 0{N n ~ 13 ) for any /3 < 1 — 2p — (k — !)/(« + 1). 
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Hence 



(|Spec (f N + N-»Q N (u)) n Q\) = / Etr(J N + N^Q N {u) - z)- l dz 

27T2 J m 

dC(w 



2 ™ ./flfl Jj2„ f[w)-z 

v ■ 1 r dz 



T 2„ 2ni J 9n f(w) - z 
N n vo\ T 2n(f- 1 (n)) + 0(N n -P 



dz + 0{N n -P) 
dC{w) + 0{N n ~P) 



which is the statement of the theorem. 



□ 



Proof of Proposition \5.1 Without loss of generality we can assume that G 7. From 



Proposition 4.1 we already know that /at (7) can be approximated by a deterministic ex- 
pression 



(5.3) 

with, if a/5N n > 



I N (-y) = / tr(/W + aiP(f N f N /a')U N V£ - z)~ L dz , 



I N {l) - In{i) = O (e- £ »»/ ri + d\og (|)) 



for some c > 0, where d is the rank of ip^f^f^/a 2 ). We choose a as in (5.1), a = h p , 
where 

1 1 

h = — - , < p < - . 
2ttN ' y 2 



In view of Proposition 2.9 d = 0(h n+ P K ") and this shows that for this choice of a and for 
S satisfying the condition in the proposition, with po > n + 1/2, 

7jv(7) - 7^(7) = 0{h~ n+e+Kp + exp(-c /i n - po+p ) = 0(|7|/r rt+(K - 1)p - £ ) . 



Thus we will prove ( |5.2 ) by showing that 



tr(f N + atfj(f N r N /a 2 )U N V^ - z)' 1 = N n / (f(w) - z)- x dC{w) 

(5.4) Jj2n 

+ 0(h- n+1 - 2p ) + 0{h~ n+p{K - X) ) . 

We first show that it is enough to consider z = 0. In fact, let Un(z)Sn(z)Vn(z)* be the 
singular value decomposition of f N — z, and put 

B N (z, w) = (f N -w + ai>((f N - z)(f N - zY/a^U^V^z))- 1 . 
Then tr (B N (z, z) - B N (0, z)) = 
a tr (B N (0, z) {ip{f n f* N / o?)U N V* N - ^{{f N - z)(f N - zf / 1 a 2 ))U N {z)V* N {z)) B N (z, z)) . 
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Since rank^((/jv - z)(f N - z)*/a 2 )) = 0(h- n+K P) for 267, and B(z,w) = Op^l/a) 
for \z — w\ < a/C, we obtain 

tr (B N (z, z) - B N {0, z)) = (D(h- n+ ^-V) , 



which can be absorbed in the error on the right hand side of (5.4). Thus we only need to 
prove (5.4) with the left hand side replaced by B(z, z) and we can simply take z = 0. 



In other words we now want to prove 
(5.5) tr(f N + ai}{f N r N /a 2 )U N V N )- 1 = N n 



dC(w) 

T 2« f(w) 



+ 0(h 



-n+l-2p\ 



+ 0(h 



-n+p(K-l)\ 



The difficulty lies in the fact that the operators fx + ^{f^f^/a^U^V^ do not seem 
to have a nice microlocal characterization. We are helped by the following identity: if 
■0 G C£°(K, [0, 1]) is equal to 1 on the support of ip then 

(1 - HF N fN/a 2 ))(f N + a^{f N r N /a 2 )U N V^)- 1 = 

(1 - 4>(r N fN/a 2 ))r N (fNf* N + a^ihf^/a 2 ))- 1 . 

This is a consequence of an identity from linear algebra: 



Lemma 5.2. Let A be a matrix and USV* be its singular value decomposition. Ifip,ip £ 
C^°(1R; [0, 1]) ; ip is equal to 1 on [—1, 1], and ip is equal to 1 on the support of ip, then 

-1 



(5.7) (1 ~^(A*A))(A + ^j(AA*)UV*) 
Proof. We first note that 



1 - iP(A*A))A*(AA* + iPiAA*))- 1 . 



and similarly vJj(AA*) 
get 

(1 - ${A*A))(A + i){AA*)UV 



A* A = VS 2 V* , 4>(A*A) = V4>(S 2 )V* , 

Uip(S 2 )U* . Since 5* is a diagonal matrix, and (1 — = 0, we 



1 = V(l - 4>(S 2 ))V*V(S + i)(S 2 )y l U* 
= V(l-4>(S 2 ))(S + ^j(S 2 ))- 1 U* 
= V(l - 4>(S 2 ))S(S 2 + tpis 2 ))-^* 

V(l - 4>(S 2 ))V*^) (VSU*) (U(S 2 + ^(s 2 ))- 1 ^) 

- ?p{A*A))A*{AA* + ip(AA*)y l , 



concluding the proof. 



□ 



The identity (5.6) follows from (5.7) by putting A = /jv/a, U = Un, and V = Vjv- Using 
this we we will find a new expression for the left hand side of (5.5) so that the identification 
with the right hand side will follow from a suitable semiclassical operator calculus. 
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Lemma 5.3. We have the following approximation for the left hand side of ( |5.5 ): 
(5.8) tr(f N + o^{f N r N /a 2 )U N V^)- x = tr r N {f N f* N + a 2 ^(f N r N /a 2 ))- 1 + 0(h,- n+ ^) . 



Proof. We use (5.6) and first note that 1 — ip can be removed from the left hand side since 



(5.9) 



tr i>(f* N f N /a 2 ){f N + a^(f N r N /a 2 )U N Vj 



* \— i 
n) 



0((rank^(/*&/^ . 
The same argument works for the right hand side once we observe that 

r N (f N r N + a^iMya 2 ))- 1 = CW(l/a) , 

and this follows from using the singular value decomposition since for non-negative diagonal 
matrices 

SNiS^ + a^iS^/a 2 ))- 1 < 1/a. 

□ 



In view of (5.5) and the lemma we have to prove 
(5.10) tr r N (f N F N + a 2 ^(f N f* N / a 2 ))- 1 = N n 1 (1C[,r) 



+ 0{h 



-H + 1-2/JN 



+ 0{h 



-n+(K—T)p\ 



, T 2n f(w) 

but that follows from the calculus developed in §|2j In fact, with the a-order function 
m(w,a) = a 2 + \f(w)\ 2 , given in Lemma 2.6, 



/n/n + « 2 ^(/iv/^/a 2 ) = T N , TeS(m,a), 
T = T Q + h l - 2 PT ll T (w) = \f(w)\ 2 + a 2 ^(\f(w)\ 2 /a 2 ), TjGSM 



where we also applied Lemma 2.8 We also have T > m/2 and hence 

1/T g S(l/m, a) , 1/Tg S(l/m, a) . 
Since / G S(y/m,a), we conclude that 

f* N (fNf* N + a^hf^/a 2 ))- 1 = P N , P G 5(l/v^, a) 

/» 



P = P + h 1 - 2p P 1 , P 1 eS(l/y/m), PoH 



|/H|* + aV(|/(«OI7a 2 )' 



We now apply Lemma 2.5 and obtain (with k^> n) 

tr f* N (fNfN + aV(/jv/jv/a 2 )) -1 = / ^H^H + 0(iV~ fc+ ") 



sup 

lfll<fe 



AT" / P Q (w)dC(w) + 0(/ i - n +( 1 - 2 ") + ^-n+fca-p)) / m(w,a)- 1/2 c/£(«;) 
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1-1 



We have m(w,a) l ^ 2 < | 1 and (1.2) at z = with k > 1 implies that \f(w)\ ' is 

integrable (k, = 1 would mean that |/(w)| _1 is in weak L 1 ): 

P /*00 /»00 

/ \f{w)\~ l dC(w)= £({\f(w)\<t})r 2 dt= 0(mm(t K ,l))t~ 2 dt<oo. 

JT 2n JO JO 

It remains to show that 
(5.11) 



\P (w) - f(wY l \dC(w) = 0(h*«- 1) ) . 



def 



Putting ip(x) = ij)(x ), we rewrite the left hand side above as 



C({\f(w)\<t})d t 
+ 



-a 2 ip(t/a) 



t(t 2 + a 2 ip(t/a 
C({\f(w)\<t}) 



dt 



c({\f(w)\<t}) Z a l it/aWit/a l di 



a 2 tp(t/a)(3t 2 + a 2 (p(t/a) + a 2 {t/a)y'{t/a)) 



t 2 {t 2 + a 2 V {t/a)Y 



t 2 (t 2 + a 2 ip(t/a)) 
dt 



< C 



2a 



t K - 2 dt = Cat 



which is (5.11). Since we have now established (5.10) this also completes the proof of 
Proposition 5.1 □ 
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